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Abstract 

We give an integral variational characterization for the speed of fronts of the nonlin- 
ear diffusion equation ut = Uxx + f{u) with /(O) = /(I) = 0, and / > in (0, 1), which 
permits, in principle, the calculation of the exact speed for arbitrary /. 
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1 Introduction 



The problem of the asymptotic speed of propagation of the interface between an unstable and 
stable state has received much attention in connection with different problems of population 
growth, chemical reactions, pattern formation and others. We refer to ^ for a recent review 
and references. The best understood of such problems is that of the nonlinear reaction 
diffusion equation 

Ut = Uxx + f{u) (la) 

with 

/(O) = /(I) = 0, /'(O) > and / > in (0, 1) (16) 

for which Aronson and Weinberger [AW] have shown that any positive sufficiently localized 
initial condition m(x, 0) evolves into a front that joins the stable state u = 1 to m = 0. The 
asymptotic speed at which the front propagates is the minimal speed c* for which there is 
a monotonic front joining m = 1 to u = 0. Moreover they show that the selected speed is 
bounded above and below by 

<c*<2 sup{\/^ I u G (0, 1)} (2) 

and that the asymptotic selected front approaches the fixed point u = exponentially with 
slope 

m = -^ic* + ^c*'-Af'm. (3) 

The lower bound cl = 2^ f'{0) is that predicted from a physical argument, the linear marginal 
stability hypothesis [||. For concave functions /, the upper and lower bounds coincide and the 
speed is exactly the linear value. However, the asymptotic speed of propagation can still be 
the linear value even when the upper and lower bounds do not coincide as explicit examples 
and a variational characterization |^ which provides improved upper bounds show. We have 
recently obtained an improved lower bound on the speed of the front that enables one to 
decide when the selected speed is greater than the linear value case which is referred to as 
nonlinear marginal stability selection. There have been several reformulations of Aronson's 
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and Weinberger's rigorous results for the nonlinear diffusion equation aiming to their heuristic 
extension to other higher order and pattern forming equations 0, 0, ^ 3- None of these 
approaches however provide the means to calculate a priori the velocity of the fronts. 

The purpose of the present work is to extend our previous result |^ to show an inte- 
gral variational characterization of the speed of the fronts of equation (1) which enables, in 
principle, its exact calculation for arbitrary /. Our main result is the following 

Theorem 1 (Variational Characterization of c*) Let f G C^(0, 1) with /(O) = /(I) = 
0, f (0) > and f{u) > for u E (0, 1). Then 



c* = max{2^/'(0),J}, 

where 

J = snp{I{g) \geV}. 

Here, 

lis) = 2&^^ (5) 

Jo gdu 

and V is the space of functions in C^{0, 1) such that g > 0, h = —g' > in (0, 1), g{l) = 0, 
and Jq g{u) du < oo. Moreover, if c* ^ 2y^/^(0), J is attained at some g eT>, and g is unique 
up to a multiplicative constant. 

In Section 2 we prove Theorem 1, and in Section 3 an example is given. 



2 Proof of the Variational Characterization 

We are interested in the calculation of the minimal speed for which equation (1) has a mono- 
tonic travelling front u{x, t) = q{z) with z = x — ct joining u = 1 to u = 0. Since the selected 
speed corresponds to that of a decreasing monotonic front, it is convenient to work in phase 
space. Calling p{q) = —dq/dz, where the minus sign is included so that p is positive, we find 
that the monotonic fronts are solutions of 

P{q) ^ - cp{q) + f{q) = 0, (7) 
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with 

p(0) = 0, p(l) = 0, p>0 in (0,1). (8) 

As shown by Aronson and Weinberger 0, the asymptotic speed of propagation of fronts of 
equation (1), c*, is the minimum value c for which there is a solution of (7) and (8). Aronson 
and Weinberger have proved that there is a (unique) p satisfying (7) and (8) for c = c* (see 
0, Section 4). Moreover, the solution p is such that p{q) ~ |m|g near q = 0, where |m| is the 
largest root of the equation 

- c*x + /'(O) = 0, (9) 



I.e., 

Iml 
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We find it convenient to introduce the parameter A defined as A = c*/|m|. In terms of A one 
can write 

It is straightforward to verify that whenever 1 < A < 2 the value of |m| given by (10) 
corresponds to the largest root of (9) and therefore to the asymptotic slope at the origin of 
the selected front [^. At A = 2 the speed c* attains the linear value c^. 
Proof of Theorem 1. The proof of the Theorem is done in two steps. First we show that 



c* > max{2v//'(0), J}, (11) 
and then we show that either c* = 2yJ /'(O) or c* = J. From the results of Aronson and 



Weinberger [AW] it follows that c* > 2^/'(0) (see equation (2) above). Thus, to prove (11) 
we need only to show that 

c* > I{9) 

for all g eT). This latter fact has been proven by us in |^. We repeat here the argument for 
completeness. Let g be any function in T>. Multiplying equation (7) hj g/p and integrating 
with respect to q we find after integrating by parts, 

^_Io{hp + ifg) dq 



Io9{(l)dq 
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(12) 



However since p, h, f and g are positive, for every fixed q 

hp+Ml>2^[f^ (13) 
P 

hence 

c* > 2 iv^ZM (14) 
/o ^c?? 

which proves (11). 

To finish the proof of the theorem notice that if c* = 2^ f'{0) we are done. Therefore, 
let's assume c* ^ 2y'/'(0) (in fact, c* > 2^f'{0)). We will show that c* = J and that there 
exists g such that c* = /(^). 

Let p{q) be the positive solution of (7) satisfying (8). The existence of such a solution has 
been established by [AW]. Moreover, p{q) ~ \m\q near ^ = 0. 

In terms of p, let us define the function v{q) as 

ri c* 

v(q) = exp{ / — dq} (15) 

■'go P 

for some fixed < go < 1- Since p is continuous and strictly positive in (0, 1), v is in (0, 1). 
Moreover, v satisfies 

u' r* 

(16) 



V p 
Now choose 

From (15) and (17) it follows that 

p{q) 

expi 

Iq p 

Clearly, ^(1) = 0, g{q) > 0, and since p e C^{0,1), g e C^(0, 1). In order to show that g eV, 
we must verify that h = —g' > in (0, 1). It follows from (17) that 

and taking the derivative of (16) with respect to q we have 



^(g) = ^exp{r-rfg}. (18) 

C* Jq V 



Mg) = /(?)exp{/ -dq}>0. (21) 



Using equations (7) and (16) to eliminate p' and v' in (19b) we find 

v" c* f 

- = -T (20) 

V p-^ 

Since v, f and p are positive, v" > 0; hence, it follows from (19a) that h — —g' > in (—0, 1). 
Alternatively, it follows from (7) and (18) that 

flo c* 
p 

Finally we must show that Jq g{q) dq < oo. As we have seen, ^ is a continuous, positive and 
decreasing function in (0, 1). Hence, it is bounded away from the origin. Thus, to determine 
whether Jq g{q) dq is finite we must study the behavior of g near q — 0. Since we know that 
p ~ \m\q near 0, it follows from (16) that 

- ^ —- - - 
V |m| q q ' 

which in turns implies that 

1 1 
Xq^ ^ 

near q — 0. Therefore, if A < 2 (i.e., if c* > cl), we have g{q) dq < oo and g eT>. 

Having verified that g E V we now show that for this choice of g the equality holds in 
(13). We must verify that 

fg 

lip -r 

or equivalently, that 



p V 



h f9 
hp — — . 

P 

This follows directly from the definition of g, in fact from (19a) hp = —g'p = v"p/v'^. 
Using equations (16), (20) and the definition (17) of g the result follows. Having shown that 
whenever A < 2 (i.e., whenever c* > cl) there exists a function ^ e X> for which we obtain 
the exact speed we have proven (4). 



Note: The uniqueness of the maximizer follows from the fact that we need to choose g in such 
a way that (13) is satisfied as an equality. Alternatively, we can prove uniqueness directly 
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from the concavity of I{g)- From the homogeneity of I{g) in g we can restrict, without loss of 
generahty, to g's such that g{q) dq = 1. Assume c* > cl and that there are two g's in P, gi 
and g2 say, such that I{gi) = /(fi'2) = c* (and J gi = / (?2 = !)• Then, let g = agi + (1 — a)fl'2, 
whith < a < 1. Clearly, J g = 1, and h = —g' = ahi + (1 — a)h2 > 0. Therefore, 

I{g) = \lf{agi + (1 - a)g2){ahi + (1 - a)/i2) 

However, by the Cauchy-Schwarz inequality. 



\agi + (1 - a)g2)^{ahi + (1 - a)/i2 > ay /^iS'i + (1 - «) V /i25'2 
with equality if and only if gi = (3hi and (72 = /5^2- Thus, 

c* > > al{g,) + (1 - a)/((72) = c* 
which implies hi/gi = h2/g2, and therefore gi = g2- 

3 Example 

In this section we illustrate the results by applying it to the exactly solvable case f{u) 
u{l — u){l + au) for which it is known that 

+ if a > 2 



2^/'(0) = 2 ifa<2 

The idea is to combine the variational characterization contained in Theorem 1, which given 
a trial function yields lower bounds for c*, with the variational characterization of Hadeler 
and Rothe which provides upper bounds on c*. 
For a given a > 0, let us take 



^(^) = A 1 A = 1 + f . 



Then 



h = -g'='^--^{\-l){l + aq) 
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and [|TT 

1 , , r(A + 2)r(2-A) _ ^ 

- 9{q)dq = ^ if A<2, 

r(4) 

i.e., if a > 2. 

We obtain 



and therefore 



V r(2-A)r(5)y 

Using the definition of A and T{z + 1) = zT{z) we obtain 
for a > 2. Therefore, from Theorem 1 we have that 



c* > 



- + J^ if a > 2 

a y 2 



2^/'(0) = 2 if a < 2 

To prove the reversed inequahty we use Hadeler and Rothe's variational characterization 

c* < sup{p(«) + 

u p[u) 

for any p G C^(0, 1) such that p(0) = and p{u) > in (0, 1). 



For a > 2 choose p{u) = u{l - u). Then p + (f/p) = y'a/2 + ^J2/a. For a < 2 

choose p = u{l —u). Hence, p+{f /p) = 2 + (a — 2)m and supo<u<i(p + (//p)) = 2. Therefore, 
from Hadeler's and Rothe's bound we get 

f \f^+\fl if a>2 
[ 2^f'{0) = 2 if a < 2 

which combined with our lower bound gives the desired result. 
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4 Conclusion 



We have given a variational characterization of the minimal speed for which the nonlinear 
diffusion equation has monotonic fronts. As [AW] have shown this is the asymptotic speed of 
propagation of a sufficiently localized positive initial condition u{x, 0). 

The variational principle we have derived here can also be used to study the dependence 
of c on the parameters of /. Monotonicity properties can be immediatly derived. Derivatives 
of c with respect to parameters of / can be obtained using the Feynman-Hellmann formula. 
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